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Recent Wilkinson Microwave Anisotropy Probe (WMAP) measurements have determined the 
baryon density of the Universe fit with a precision of about 4%. With fit tightly constrained, 
comparisons of Big Bang Nucleosynthesis (BBN) abundance predictions to primordial abundance 
observations can be made and used to test BBN models and/or to further constrain abundances 
of isotopes with weak observational limits. To push the limits and improve constraints on BBN 
models, uncertainties in key nuclear reaction rates must be minimized. To this end, we made 
new precise measurements of the 2 H(d,p) 3 H and 2 H(rf, n) 3 He total cross sections at lab energies 
from 110 keV to 650 keV. A complete fit was performed in energy and angle to both angular 
distribution and normalization data for both reactions simultaneously. By including parameters 
for experimental variables in the fit, error correlations between detectors, reactions, and reaction 
energies were accurately tabulated by computational methods. With uncertainties around 2% ± 1% 
scale error, these new measurements significantly improve on the existing data set. At relevant 
temperatures, using the data of the present work, both reaction rates are found to be about 7% 
higher than those in the widely used Nuclear Astrophysics Compilation of Reaction Rates (NACRE) . 
These data will thus lead not only to reduced uncertainties, but also to modifications in the BBN 
abundance predictions. 

PACS numbers: 26.35.+C 25.10.+S 25.45.-z 25.60.Pj 



I. INTRODUCTION 

The standard Big Bang model explains remarkably 
well many features of the Universe which are otherwise 
difficult to reconcile. The standard Big-Bang nucleosyn- 
thesis (BBN) model consists of a small network of nuclear 
reactions occuring at energies easily obtained in the lab. 
The outcome of Standard BBN is determined almost en- 
tirely by the nuclear reaction rates and the baryon to 
photon ratio, 77, of the Universe, which is directly related 
to the baryon density fib {flbh 2 = 3.66 x 10 7 ?7, where h is 
the Hubble constant in units of 100km xs _1 x Mpc^ 1 Q ). 
With detailed network calculations and knowledge of the 
nuclear cross-sections, all of the primordial abundances 
can be precisely calculated as a function of 77. Knowl- 
edge of any one abundance or a separate determination 
of 77 thus allows all other primordial abundances to be in- 
ferred from the Standard BBN model. Knowledge of any 
two of these observables produces a check of the model 
itself. 

In the past, quantitative understanding of BBN has 
been limited by uncertainties in the observed primordial 
abundances and the value of r\. Until recently 7/ was 
treated as a free parameter. As observations have im- 
proved, the value of r\ has become a well determined in- 
put and is instead used, along with the nuclear reaction 
rates, to predict the primordial abundances. Regardless 
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of which values are assumed and which are predicted, the 
uncertainties of the cross sections themselves are becom- 
ing a significant factor in precision tests of BBN models. 
Improved analyses depend on having accurate and precise 
knowledge of the nuclear reaction rates, their uncertain- 
ties, and the correlations in those uncertainties. 



A. Primordial Deuterium Observations 

Recent measurements of absorption lines in high red- 
shift, metal poor, QSO-back-lit gas clouds have con- 
strained the primordial deuterium abundance D to the 
impressive interval of D/H = 2.78±^ x 10" 5 0, ex- 
pressed relative to the hydrogen abundance H. The mea- 
surement and analysis procedure is well described in the 
review article of Tytler et al. Q. Currently systematic 
scatter limits the precision of the deuterium abundance 
observations 0, but as more data arrive and the sys- 
tcmatics become better understood this could quickly 
change. 

Measurements of primordial deuterium abundances 
0, , along with recent measurements of r\ from WMAP 
data 0, , bring the nuclear reaction rates increasingly 
closer to being the limiting factors in testing the consis- 
tency of the Standard BBN model. 



B. Cosmology Enters the Lab 

The value of r\ determined from the WMAP results 
is in good agreement with primordial deuterium abun- 
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dance measurements of Kirkman et al. 2] . However lim- 
ited amounts of data and the systematic uncertainties in 
abundance measurements and in reaction rates have lim- 
ited the level of precision at which the BBN models can 
be tested to around 10% for most observables. 

With the uncertainties of previously existing data, the 
2 H(d,p) 3 H and 2 H(d, n) 3 He reaction cross-sections at 
energies in the range of a few hundred keV make large 
contributions to the uncertainties in deuterium abun- 
dances as predicted by network calculations |sl. Ijj. [ToL fTl| ■ 
Nollett and Buries |8| provide sensitivity functions es- 
timating contributions to the deuterium and 7 Li abun- 
dances of several reactions as a function of energy. The 
most relevant energy range of both the 2 H(d, n) 3 Re and 
2 H(d, p) 3 H reactions extends from roughly Ed — 100 keV 
up to about Ed = 700 keV. Prior to our measurements, 
high precision data for these reactions were very limited 
in this range, as shown in Fig.^ The 2 H(<i, p) 3 H reaction 
data is qualitatively very similar. From Ed — 325 keV 
to the top of the BBN energy range there were very few 
data points, all having large uncertainties. Even at lower 
energies of significance to BBN, the uncertainties of pre- 
viously existing data were around the 10% level or only 
slightly better. For the energy range above 325 keV, one 
of the more complete data sets was that of Ganeev et 
al. |l2| . These data are not included in the NACRE data 
compilation considered the most prominent collec- 
tion of experimental rates for reactions of astrophysical 
significance. 
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FIG. 1: (Color) Prominant data sets for 2 H(d, n) 3 He cross- 
sections at BBN energies 

In order to aid in tightening the BBN constraints, 
we measured total cross sections for both of these deu- 
terium burning reactions at lab energies ranging from 
about 112 keV to 646 keV. By carefully establishing ex- 
perimental procedures in order to maximally cancel the 
dependencies of yields on experimental parameters, we 
have obtained about 2% statistical uncertainties plus or 
minus a 1% systematic scale uncertainty. Furthermore 
our \ 2 analysis techniques have allowed optimal use of 



the data and provided correlations in the uncertainties 
of the two reactions across the range of energies. The 
new data form a significant improvement of the inputs to 
the BBN calculations and facilitate an emerging era of 
high-precision BBN. 



II. EXPERIMENTS 

To determine the integrated cross sections for the 
2 H(<i,p) 3 H and 2 H(d, n) 3 He reactions, relative angular 
distributions of the differential cross sections were 
measured, henceforth denoted as simply a{0), or a(E : 6), 
where E and 9 represent the energy and reaction an- 
gle respectively. We normalized these distributions by 
measuring absolute differential cross sections at selected 
fixed angles. This general procedure was performed for 
eight energies, £^=120, 180, 240, 320, 390, 480, 560, and 
650 keV. These are approximate nominal energies and 
are the values used here to refer to the various data sets. 
The determination of the exact energies will be discussed 
in Sec. Ill Cl Generally data for both reactions were ob- 
tained simultaneously. 

All data were taken using the TUNL Low-Energy 
Beam Facility (LEBF) and the High- Voltage Target 
Chamber The LEBF is composed of the Atomic 

Beam Polarized Ion Source [l£j and Mini- Tandem ac- 
celerator The combined acceleration potentials 
of the source, Mini- Tandem and HV chamber provide 
deuteron and proton beam energies up to approximately 
680 keV. All targets used were carbon-based, self- 
supporting transmission targets. 



A. Normalization Technique 

For the present measurements, the central technique 
used to determine the absolute cross sections was to com- 
pare them to those of a reference reaction. We used p-d 
elastic scattering for which multiple absolute measure- 
ments of differential cross sections exist at ener gies near 
the range of interest for our measurements (l7l JlH. . 
The data of Brune et al. was taken with the explicit 
purpose of normalizing relative differential cross-section 
data taken with the same equipment as the present work. 
All of these measurements agree, to within about 1%, 
with theoretical few-body calculations with no free pa- 
rameters ^3 US • We used the calculated values as the 
cross-section reference for our measurements. 

In order to compare yields for the d-d interactions 
and p-d interactions under similar conditions, we alter- 
nated between deuteron beam and proton beam inci- 
dent on the same self-supporting, deuterated carbon tar- 
gets. The targets were produced by plasma deposition of 
fully-deuterated methane gas and were typically about 
30 /ig/cm 2 thick, with roughly equal numbers of carbon 
and deuterium atoms. 
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The ratio of cross sections for the two observed inter- 
actions is given by 

a d (E d ,9) _ riptpAQp N d 
a p (E p ,9) n d t d Att d N p ' 

where n is the incident number of beam particles, N is the 
detected number of particles, and t is the areal density of 
target nuclei. The subscripts p and d refer to measure- 
ments with proton and deuteron beam respectively. By 
using the same experimental setup and target for both 
beams, the solid angles and target thicknesses were made 
to cancel. 

One of the larger obstacles in a direct cross-section 
measurement is measuring the time- integrated beam flux 
incident on the target. For the energies of the present ex- 
periment, the determination of the number of incident 
beam particles by charge collection and integration is 
made difficult by angular straggling and charge exchange 
in the target. 

In order to determine the ratio of proton beam parti- 
cles to deuteron beam particles, a thin layer of gold, ap- 
proximately 10 or 20 A thick (lA of gold is equivalent to 
1.7 x 10 15 atoms/cm 2 ), was evaporated onto one surface 
of the target. The gold layer was then oriented on the 
upstream side of the target, facing the incident beam. 
A back-angle detector was used to monitor Rutherford 
elastic back-scattering from gold for both proton and 
deuteron beams. 

The ratio of the number of incident beam particles for 
the two beams is given by 

nP = N* <r Au{dtd) {E d ,0) 

Nda^p^iEp^)' U 

where OAu(d,d) (E d , 0) and ctau(p, p ) (E p , 0) are respectively 
the differential cross sections for proton and deuteron 
scattering from gold at the same lab angle. 

When beam is incident on the amorphous deuterium 
targets, deuterium is depleted from them, albeit slowly. 
If a target is depleted in a spatially non-uniform manner 
then two beams striking slightly different areas at differ- 
ent times may not interact with a target of the same 
thickness. High intensity beam currents tend also to 
produce macroscopic defects in targets such as rips and 
holes, especially at these energies where energy loss in the 
target is large. Depletion and stability problems seem to 
depend not only on the total number of particles which 
pass through the target, but also on the rate at which 
they pass. 

A constant or slowly varying target thickness must be 
maintained between proton and deuteron runs. Since 
the reactions of interest have cross sections on the or- 
der of mb/sr, as opposed to b/sr for clastic scattering, a 
much larger integrated beam flux was required in order 
to acquire sufficient reaction data. Such high beam fluxes 
would produce prohibitively large target variations. 



To avoid this problem, we first normalized d-d elastic- 
scattering yields to the p-d yields using small beam cur- 
rents, below 10 nA. At these currents target depletion 
was almost imperceptible, no more than 1% per hour, and 
many targets were structurally stable for several hours. 
Data was taken later with higher beam currents, around 
50 to 100 nA, and the ratio of d-d elastic-scattering yields 
to 2 H(d, p) 3 H and to 2 H(rf, n) 3 He reaction yields was ob- 
tained. Since the elastically scattered deuterons and the 
reaction products were measured simultaneously in this 
last step, they passed through the same target thick- 
nesses. This final step adds negligibly to the overall sys- 
tematic error of the measurement. To further ameliorate 
the depletion problem, the beam switching technique was 
employed several times over short intervals, and small 
entrance collimators were used to define precisely the in- 
cident beam trajectory. 

In order to change beam types rapidly, we injected a 
mixture of deuterium and hydrogen gas into the ECR ion- 
izer of the Polarized Ion Source [TJj . With this dual-beam 
source configuration, by changing only the inflection- 
magnet current, we could easily put over 100 nA of either 
beam on target. 



B. Experimental Setups and Procedures 

The HV chamber which was used for these measure- 
ments is described in Ref. 01- The chamber's basic fea- 
tures are illustrated in Fig. [3 Two independently rotat- 
ing plates are installed, one on the top of the chamber 
and one on the bottom, which can be used to mount 
detectors on the left and right side of the chamber re- 
spectively. The fixed monitor detectors are placed above 
and below the reaction plane having a view of the target 
which was not obstructed by the rotating detectors. The 
beam enters the chamber through an acceleration tube 
allowing for an acceleration through a 200 kV potential. 
It then passes through vertical and horizontal entrance 
slits which define the beam position. 



1. High-Energy Normalization Setup 

For the highest four beam energies, 390, 480, 560, and 
650 keV, wc used the procedure described in Sec. Ill Al 
to normalize the differential cross sections. Seven high- 
resolution, ion-implanted silicon detectors were posi- 
tioned in the chamber to observe elastic scattering from 
deuterium. All detectors had either 300 or 500 /im de- 
pletion depths, and 1.50 cm 2 active areas. Four were 
placed on the rotating detector tracks and set at 35.0° 
and 48.0° on both the left and right sides. Two more were 
mounted as out-of-plane monitors at 27.5°. One detec- 
tor was placed on the top rotating plate at 125.0° left 
for the purpose of monitoring Rutherford backscattering 
from gold as explained in Sec. Ill Al The rotating detec- 
tors were placed approximately 18 cm from the target 
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FIG. 2: (Color online) Schematic of the experimental set-up. 
The entire scattering-chamber can be raised to a potential of 
±200 kV to accelerate or decelerate the beam. 
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FIG. 3: Spectrum from proton elastic scatterring from a 
deuterated target at f9i ab = 27°, E p =560 keV. 
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FIG. 4: Spectrum of d+d reaction products at #i a b = 35°, 
E d = 480 keV. 



with collimators directly in front of them with openings 
0.46 cm wide and 0.95 cm tall. The out-of-plane and 
back-angle monitor detectors were mounted at approx- 
imately 30 cm and 10 cm from the target respectively. 
The out-of-plane monitors were given collimators 0.64 cm 
wide by 0.95 cm tall, and the back-angle monitor was 
equipped with a circular collimator of 1.27 cm diameter. 
Cylindrical aluminum tubes were placed in front of all 
detector collimators to help reduce possible background 
scattering from various sources such as the entrance slits. 
Small permanent magnets were placed around the tubes 
to reduce the effects of electrons released from the target. 
Systematic cancellations inherent in the procedure re- 
moved the need for precise knowledge of collimator sizes 
and radii. An aluminum plate serving as a Faraday cup, 
was placed at the back of the chamber to integrate the 
beam current. Current integration was used for beam 
adjustment and diagnostic purposes only. 

The two pairs of detectors on the rotating plates were 
used to observe p-d elastic-scattering yields during pro- 
ton beam runs. A spectrum is shown in Fig. [3J The 
out-of-plane monitors were used to measure d-d elastic- 
scattering yields during the interleaved deuteron runs. 
The in-plane detectors were again used for observing re- 
action products (see Fig. 0} during the reaction runs 
while simultaneously measuring d-d elastic-scattering 
yields in the out-of-plane monitors. Using different de- 
tectors for the two simultaneous deuterium-beam mea- 
surements provided more flexibility for the choice of an- 
gles used and more easily facilitated differing signal gain 
requirements for the elastic-scattering and reaction prod- 
ucts. Solid angle cancellations required that all d-d elas- 
tic yields were measured in the same dedicated detectors 
and that the reaction yields were measured in the same 
detectors as the p-d elastic scattering yields. 

Measurements of the d-d scattering yields were com- 
plicated by the presence of hydrogen contamination in 
the targets. Within the available increments, the an- 
gle of 27.5° was experimentally determined to produce 
the best spectral separation of the d-d scattering peak 
from both the 1 H(d,p) and 12 C(d, d) peaks. At £i a b = 
390 keV, the lowest energy where the beam-switching 
normalization was performed, adequate spectral separa- 
tion of the d-d scattering peak from the 1 H(d,p) peak 
could not be achieved. Instead the yield from the well- 
separated 1 H(rf, ri) scattering peak was used to determine 
the 1 H(ri,p) yield which could then be subtracted from 
total yield of the unresolved 2 H(rf, d) and 1 H(d, p) peaks. 
If the ratio of hydrogen to deuterium were constant in 
every target, the 1 H(d,p) yield would introduce no er- 
ror since it would still provide a consistent, proportional 
measure of the amount of deuterium in the target. Typi- 
cally the 1 H(ri, p) yield was less than 10% of the 2 H(d, d) 
yield. The ratio of the two varied by about 30% from tar- 
get to target, thus producing about a 3% error. Only a 
rough correction was needed to reduce this to a negligible 
effect. In order to measure the ratio of the 1 H(<i, d) yield 
to 1 H(d,p) yield, monitor spectra were collected at sev- 
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eral energies for deuteron scattering from a hydrogenated 
carbon target. 



2. Low-energy Normalizations 

It was not possible to normalize data at the lowest 
energies, 120, 180, 240 , and 320 k eV using the method 
described in Sees. Ill Al and III B II At these energies it 
was not possible to resolve all needed elastic peaks with 
sufficient precision. We instead normalized reactions at 
these energies directly to reactions at 480 keV using the 
differential cross sections already obtained for 480 keV 
from the previous method. 

For this normalization, four detectors on each side of 
the chamber were placed at 60°, 44°, 28°, and 13°. As be- 
fore, one detector was placed at 164° on the top rotating 
plate to monitor backscattering from an upstream gold 
target layer. All detectors were placed approximately 
11 cm from the target and collimated with 1.27 cm di- 
ameter circular collimators. 

At these low bombarding energies, amorphous deuter- 
ated carbon targets usually became very fragile, and were 
replaced with deuterated parapolyphenol (DPP, chemical 
formula: C6D4 ) targets [21j . We constructed the targets 
by evaporating an approximately 5 /ig/cm 2 -thick layer of 
DPP onto a 5 ^g/cm 2 carbon foil. As in the high-energy 
normalizations, a roughly 3 /xg/cm 2 -thick Au layer was 
evaporated onto the DPP surface. Again the gold layer 
faced upstream so that the beam first passed through 
the gold, then the DPP and then the remaining carbon 
backing. 

The actual normalization procedure was fairly 
straightforward. We bombarded the target with deu- 
terium beam at 480 keV and observed ratios of reaction 
yields to Au backscattering yields. We then changed the 
beam energy to the energy of interest while using the 
same target. The procedure was repeated several times. 
Taking the ratios of the two resulting sets of normalized 
reaction yields and dividing out the two Rutherford cross 
sections provided the ratios of differential cross sections 
at 480 keV to the cross sections at each of the lower en- 
ergies. 



3. Angular Distributions 

The normalization procedures described above were 
collectively the most challenging part of the cross section 
measurements, but they only determined the differential 
cross section at a few fixed angles. To measure the rel- 
ative angular distribution of differential cross sections, 
the chamber was set up with six detectors placed on the 
rotating tracks, three on each side, with 13° separation 
between each detector. Two pairs of out-of-plane moni- 
tor detectors were placed at approximately 10° and 40°. 
Detectors were placed roughly 18 cm from the target with 



circular collimators of 1.27 cm in diameter. Beam cur- 
rents of 50 to 100 nA were used. Amorphous carbon tar- 
gets were used at the highest energies for their high deu- 
terium content. At lower energies carbon backed DPP 
targets were used for their superior durability. Using this 
arrangement, data at all angles could be normalized to 
data from the fixed monitor detectors, thus dividing out 
all beam current and target thickness information and 
leaving only a relative angular distribution of a (6) . 

The elastic-scattering rate from carbon present in the 
target was much higher than the rate of d + d reactions. 
With the high Q-values of the reactions, 4.03 MeV and 
3.27 MeV for 2 H(d,p) 3 H and 2 H(d, n) 3 He respectively, 
it was possible to separate the reaction peaks from the 
lower-energy elastic-scattering peaks for most reaction 
angles and beam energies. To avoid overloading the data 
acquisition system, Mylar stopping foils from 1 to 6 /im- 
thick, depending on the angle and energy, were placed in 
front of the detectors. 



4- Detector Electronics 

As described in Ref. 0], the High- Voltage Chamber 
accommodated many silicon surface-barrier detectors in 
various moveable configurations. For many parts of the 
experiment, the chamber was electrically isolated from 
ground and brought to a high negative potential of as 
much as 200 kV. This required that all energy signals be 
sent from the chamber via analog fiber-optic transmit- 
ters as described in Ref. j2^. Energy signals from the 
fiber-optic receivers were then routed into six Northern 
ADC's. Pulses of known rate were sent through the same 
electronics in order to measure the deadtimc. 



C. Energetics 

To achieve the desired uncertainties in the cross sec- 
tions, reaction and scattering energies must be known 
very well. The incident beam energy was determined by 
the potentials on the ion source, on the Mini Tandem, 
and on the High- Voltage Chamber. The energy calibra- 
tion of the LEBF system was previously determined by 
a procedure described elsewhere Small corrections 

were made for energy losses in the Mini Tandem carbon 
stripping foil and for the potential on the cesium in the 
charge exchange canal in the ion source. The uncertainty 
in the incident beam energy was less than 1 keV. 

Energy losses in the targets must also be well known 
in order to determine precisely the reaction and scat- 
tering energies. In particular the elastic-scattering cross 
sections used to determine the normalizations are very 
sensitive to energy, especially at the lowest energies. For 
these low-energy normalizations a three-layer target was 
used. The deuteron beam first passed through a thin 
layer of gold followed by a layer of DPP and finally passed 
through a carbon backing. As explained in previous sec- 
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tions, elastic-scattering yields were measured from the 
gold and the carbon, and reactions were measured from 
the deuterium in the DPP layer. Since energy was lost 
continuously throughout the thickness of the target, the 
reaction energy was less than the gold scattering energy 
which was less than the incident beam energy. Both in- 
teraction energies need to be known. 

By rotating the multilayer-targets by 180° and ob- 
serving shifts in the energies of elastically backscattered 
beam particles, as well as by monitoring relative scatter- 
ing yields of various targets, it was possible to measure 
and continuously monitor target thicknesses and energy 
losses. The energy losses, due primarily to the carbon 
content of the targets, did not change significantly over 
time. Final values of incident beam energies and total en- 
ergy losses in the gold and deuterated layers are given in 
Table HJ The high energy normalizations were performed 
on two-layer targets of gold and a deuterated amorphous 
carbon. The energy losses in these targets are shown in 
Table HU The impact of the energy uncertainty is dis- 
cussed in Sec. IHI E 31 



III. DATA ANALYSIS 

The goal in analyzing the data was to make full use 
of all the interconnected data sets constraining various 
physical and experimental parameters without double 
counting any statistics, and to understand the correla- 
tions in the uncertainties in the final results. By using a 
X 2 analysis to fit all the data to a model constructed from 
the relevant parameters, it was possible to maximize the 
amount of information used while appropriately handling 
and quantifying correlated information and uncertainties. 

Before performing the global fit of the complete data 
set, some initial analysis was performed on the data from 
each storage sequence, referred to as a run. The energy 
spectra were first analyzed to determine the yields in all 
peaks of interest. When required, fits were made to Gaus- 
sian peaks summed with linear or quadratic backgrounds. 
Data were then adjusted for deadtime corrections, mon- 
itor normalizations, and cross-section normalizations as 
needed. In this way a preliminary conventional analysis 
could be performed, and parameters which changed with 
each run, primarily including target layer thicknesses and 
the number of incident beam particles, could be removed 
from the data sets. 



A. High-Energy Normalizations 



simplicity we will assume that each yield was observed in 
only one detector and one data run. The detector used to 
measure p-d elastic-scattering yields will be referred to as 
detector A. The d-d elastic-scattering yield was observed 
in detector B in a beam switching run. These two runs 
used the same target, which we will call target 1, having 
a deuterium target thickness t d \ and a gold target thick- 
ness t g \. Finally d-d elastic scattering was measured in 
detector B simultaneously with the d-d reactions in de- 
tector A. These measurements were made using a differ- 
ent target, target 2. For all three measurements, elastic 
scattering on gold was measured in detector C. The su- 
perscript p indicates a proton-beam run, d a deuteron 
beam-switching run, and finally r denotes reaction runs 
with deuteron beam. These label the yields, the beam 
energies, and the integrated beam flux on target. Sub- 
scripts indicate reactions where needed. The yields of 
interest arc as follows: 



N *H(p,p) 
Au(p,p) 
N '^H(d,d) 



2 H(<2,<2) 



nPt gl An c a Auip ^(EP,6 c 
n d t dl An B <T2 H{did) (E d ,6 B 

n d t gl An C (JAu(d,d)(E d ,6c 
ri r t d2 AQ A (T2- ti{d ^(E r , 9 A 
n r t g2 AVL B (KYi(d,d) (E r , B 



(3) 
(4) 
(5) 
(6) 
(7) 
(8) 



The solution for the 2 H(d,p) differential cross section 
is given by the following: 



c 2 -R(d,p)(E r ,9 a) = o-2 H ( pp) (iJ p , 0)x 

N m(d, P ) N m(d,d) N lu(p,p) VAu(d,d){E d ,9 c ) 
N m(d, d) NL {d , d) N m (p , p) ^Au {p , p) (EP,9c) 



(9) 



In order to analyze the data one run at a time, we break 
this solution into the normalization factors a correspond- 
ing to the three different run types: 



a' 



a 



N ^K( P , P ) (J Au(p,p){E p ,9 C ) 



N d 
N d 

iV2 H(d,p) 



°Au(d,d) 



(E d ,9, 



C 



Nl 



(10) 

(11) 
(12) 



H(d,d) 



As explained in Sec. Ill Al the highest-energy differ- 
ential cross sections were normalized to the known p-d 
elastic-scattering cross sections. Here we outline the pro- 
cedure needed to obtain absolute reaction cross section 
from the p-d cross sections. 

First we construct the specific yield equations for all 
reactions and run types required for the calculation. For 



Now we have simply 

a* H{d , p) (E r ,e A ) = ^. (13) 

By substituting the yields with the systematic param- 
eters on which they depend, these factors can also be 
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TABLE I: Energy corrections for low energy normalizations: The final interaction energy is calculated using the total energy 
loss in gold and half of the DPP energy loss. The values shown here for 480 keV correspond to the target used to normalize 

the 120 keV data to the 480 keV data. 

Nominal Incident Deuteron Total Loss in Total Loss in Central Deuterium- 
Energy (keV) Energy (keV) Gold (keV) DPP (keV) Layer Energy (keV) 



480 477.0 0.5 7.1 473.0 

320 317.3 0.2 5.1 314.6 

240 237.2 0.5 7.9 232.8 

180 177.3 0.5 8.1 172.8 

120 116.8 0.4 8.5 112.2 



TABLE II: Energy corrections for high-energy normalizations: Each row describes deuteron beam energies of the data being 
normalized and the proton energy used to normalize it. Proton beams having two different energies were used on the same 
target to normalize the 560 keV deuteron data. Variations in energy losses between rows come from varying stopping powers 
and variations in target thicknesses. 

Nominal Incident Deuteron Central Deuteron Incident Proton Central Proton 
Deuteron Energy Beam Energy Beam Energy Beam Energy Beam Energy 



(keV) (keV) (keV) (keV) (keV) 

~650 653^2 646T 65315 6491T 

560 564.4 557.3 653.6 649.4 

560 564.8 560.2 

480 477.0 470.2 564.8 560.8 

390 387.4 379.2 564.8 560.4 



written in the following way: 



H(<M)( ' 



(14) 
(15) 
(16) 



which makes evident the solid angle and target thickness 
cancellations in Eq. 

To calculate a p , differential cross-section values were 
needed for p-d elastic scattering. Theoretical calcula- 
tions were provided by Kievsky et al. [2(j [23| for several 
selected energies in the range of interest for the present 
work. Two-dimensional polynomial interpolations were 
used to obtain values for the precise angles and energies 
required. The error associated with the interpolation is 
negligible. 

All cross sections for scattering on gold used in our 
analysis were calculated from the Rutherford scatter- 
ing formula with electron-screening corrections obtained 
from Ref. [24|. At the lower energies, the screening cor- 
rections were about 1% of the the values. 

The error due to target thickness variation can be seen 
by plotting the values of a p and a d as a function of run 
number. Run numbers were incremented consecutively 
as beams were switched between proton and deuteron 
beams. The results for the 650 keV normalization are 
shown in Fig. [S] The two factors, having different phys- 
ical meanings, are shown in relative units. The relative 
stabilities of the two measurements is apparent. 
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FIG. 5: (Color) Relative yields measured in beam-switching 
normalization runs: a p and a d are the normalized 2 H(p,p) 
and 2 H(d, d) yields as defined in Eqs. HOI and llll 



B. Angular Distributions 

Deadtime-corrected yields for all rotating detectors 
were normalized by dividing them by an appropriate 
combination monitor-detector yields. The resulting nor- 
malized data was sensitive only to the differential cross 
sections of the reactions and an overall normalization fac- 
tor. The relative solid angle normalizations of the rotat- 
ing detectors could be determined from cross-calibration 
data taken for that purpose, but ultimately these nor- 
malizations were left as free parameters as discussed in 
Sec. IIII Dl With four monitor detectors there was suffi- 
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cient redundancy to use an automated algorithm to re- 
ject monitors which were blocked or partially blocked in a 
particular chamber geometry. Angle and solid-angle dif- 
ferences between the four monitors required the relative 
normalizations of the monitor yields to be determined 
for every energy. The full data set of all reaction runs at 
each energy was used for this determination. Iterations 
were made alternating between the relative normalization 
procedure and the monitor rejection algorithm. 

Fig. shows the resulting angular distribution for the 
2 H(d,p) 3 H reaction at E d = 650 keV with a fit to even- 
order Legendre polynomials in center-of-mass angle. 
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FIG. 6: (Color online) Normalized angular distribution of 
2 H(rf,p) 3 H differential cross section at Ed = 650 keV. The 
data points shown are the weighted means of all data runs 
at the same angle. The curve is a fit of even order Legendre 
polynomials up to P4. 



C. Low Energy Normalizations 

Once normalized angular distributions were derived at 
E d = 480 keV, the energy cross-normalization data could 
be used to normalize the low energy angular distributions 
at fixed angles. Again, for simplicity, we will consider 
just one reaction detector R. Only the 2 H(rf,p) 3 H and 
2 H(d, n) 3 He reaction peaks were useable for the cross- 
normalizations. The reaction yields in R were normal- 
ized only to the elastic-scattering yields on gold in the 
Rutherford-scattering monitor, again labeled detector C, 
and to the appropriate Rutherford cross section. The 
normalized yield a x is defined as 

N d 

aX= ^^ aAu(dd){E \e c ). (17) 

iV Au(d,d) 

This ratio can also be written as: 

a x — ^^pcr2 H ( di p)3 H (-E d , 0b) , (18) 



using the same notation as in previous sections. 

If this ratio is determined for data taken at two 
deuteron energies, 480 keV, and for example, 120 keV, 
using the same target, then we find 



a x (E d = 120 keV) = a ddp (E d = 120keV,c?) 
a* (E d = 480 keV) <r ddp (E d = 480 keV, (?) 



(19) 



Low-energy 2 rI(d, 3 He)n reaction cross sections can 
then be determined by relating the 2 H(d, 3 He)n and 
2 H(d,p) 3 rI yields in the plentiful angular distribution 
data sets. 



D. Integrated Cross-Section Global Analysis 

The framework described in Sec. Illll is sufficient for de- 
termining total-cross section results for both reactions at 
all observed energies. However, in the previous sections 
the situation was oversimplified by considering that each 
data type was observed using one detector and consider- 
ing only one reaction at a time. In this algebraic frame- 
work and without these simplifications, it is very difficult 
to fully utilize all the data from all detectors in the nor- 
malization process while accounting for the correlations 
involved. 

To list one example, p—d elastic-scattering data and the 
normalization reaction data exist for two detector sets at 
two different angles. In order for solid angles to vanish in 
the formalism of Sec. II I II the same detector set must be 
used in the analysis of the reaction data as for the proton 
elastic-scattering data. Thus data from each detector set 
must be treated as an independent measurement. This 
actually produces complications. The fit to the angu- 
lar distribution data does not directly produce discrete 
values as a function of angle; it produces Legendre co- 
efficients and errors in those Legendre coefficients. The 
values derived for the differential cross sections at par- 
ticular angles have uncertainties which are highly corre- 
lated through the Legendre parameterization. Such cor- 
relations in the analysis inputs ultimately complicate the 
calculation of the uncertainty of a single result. 

We were also interested in calculating multiple related 
results, cross-sections for different reactions and differ- 
ent energies. The results themselves have uncertainties 
which are correlated with each other. This complicates 
the use of these values in future calculations such as the 
deuterium abundance prediction of a BBN network calcu- 
lation. In the present case, strong correlations in energy 
exist because the low energy data were normalized di- 
rectly to high energy data. Strong correlations between 
cross-section uncertainties for the two different reactions 
exist because these different reaction yields were normal- 
ized to the same data. The uncertainty correlations in 
the cross sections must be quantified. 

We have addressed these correlation issues by avoiding 
an algebraic analysis as much as reasonably possible. We 
performed a % 2 fit of the energy-dependent differential 
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cross sections. These were fit to all the normalization 
and angular distribution data for both reactions at all 
energies, all in a single and simultaneous fit. This pro- 
cedure determined uncertainties for all parameters using 
all available data which may directly or indirectly affect 
their values, ft also allowed the correlations in the un- 
certainties in those parameters to be measured using the 
usual parameter variation techniques of \ 2 fitting. 

To perform this fit we produced a realistic parameter- 
ization to describe all segments of the experiment. The 
normalized data values which were input into the fit were 
the a's of Eas. If OlO and IT7I The fit parameters repre- 
sented the unknown values in the physical representa- 
tions given in Eqs. 1141161 and El The fit thus included 
a parameterization of the energy dependent differential 
cross sections along with all needed systematic param- 
eters specific to the given data types. Actually, many 
of these variables which cancelled out of the algebraic 
analysis, such as target thicknesses, were individually un- 
determined. In most cases the appropriate fit parameters 
were carefully chosen to represent well determined prod- 
ucts and ratios of these variables. We will define here 
only the parameters of interest. The 2 H(d, n) 3 He and 
2 H(rf,p) 3 H differential cross sections can be described in 
the center-of-mass frame in terms of even Legendre poly- 
nomials P n (9) with coefficients a n ^k for the n'th Leg- 
endre polynomial, the i 'th energy and the k 'th reaction. 
To improve the minimization behavior we defined the pa- 
rameters b n ^k such that 

b n , ik = ^. (20) 
ao,ik 

The differential cross sections are then parameterized as 
<Jk{Ei,0 c . m .) = 

M (21) 

m=l 

where 2M is the highest order of Legendre polynomial 
used. The integrated cross sections are then given by 

o- fc tE i )=47ra ,«. (22) 

We also performed a fit to a parameterization which 
is continuous in energy. This has many advantages, but 
also some complications. We will not present the full 
results or description of that fit here; it is described in 
detail in Ref. Hj. 

E. Error Analysis 

The uncertainties for the normalized data input into 
the fits were calculated by first-order propagation of the 
uncertainties in the peak sums. Correlations arising from 
shared monitor normalization counts were not taken into 
account because the monitor yields generally contributed 



a relatively small amount to the overall statistical uncer- 
tainties. The final uncertainties and correlations in the 
cross-section results are determined directly from the er- 
ror matrix. We will explain this here in more detail and 
will discuss procedures applied to quantify data scatter 
and certain systematic errors. 

1. The Error Matrix 

The total cross section at each energy corresponds di- 
rectly to one of the fit parameters. The la uncertainty 
q in the cross section is simply the uncertainty in the 
corresponding parameter multiplied by a factor of 47r: 

= 47r?(oo,«) , (23) 

where T(do,ifc) is given by the square root of the diagonal 
element of the covariant error matrix corresponding to 
the parameter ao,jfc. 

2. Quantifying Scatter 

If for some fit, \ 2 P er degree of freedom, xti is n °t 1; 
then error matrix has little meaning. Some amount of 
unaccounted-for scatter is expected in the data for vari- 
ous reasons, the largest coming from fluctuations in the 
target thicknesses during the beam-switching procedure. 
In the end, xl f° r the entire fit is just over 2. 

We implemented a procedure to address this issue by 
quantifying the scatter in the data and adding an appro- 
priate amount of uncertainty to the data. This was done 
by iteratively adding in quadrature enough uncertainty 
to points in each of the data sets in order to make \ 2 P er 
datum equal to 1 for that data set. By adding a constant 
fractional error rather than multiplying the uncertain- 
ties, we were able to make better use of the full data set. 
This is a result of not over de-emphasizing points hav- 
ing uncertainties which were already large compared to 
the missing uncertainty. For a single parameter measure- 
ment, in the case where the fraction of the uncertainty 
which was initially quantified approaches zero, this addi- 
tive method reduces to the familiar technique of measur- 
ing the standard error in the mean. 

3. Angle and Energy Uncertainties 

The effects of detector-angle and energy uncertainties 
were studied by perturbing the input angle and energy 
values of all data in various ways and then reproducing 
the entire analysis for each set of perturbations. This is 
similar in principle to the uncertainty analysis performed 
by the fitting software, but for technical reasons these 
parameters could not easily be fully parameterized within 
our analysis framework. 

The detector angles were calibrated and set with an 
uncertainty of 0.1°. Fractional uncertainties arising from 
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angle determinations were found to be on the scale of 
2 x 10 -3 or less and were neglected. Uncertainties arising 
from energy determinations are larger, as much as one or 
two percent at the lowest energies, and fall to a small 
fraction of a percent at the highest energies. 

In order to consolidate these uncertainties into the er- 
ror matrix, the effect of the energy perturbations on the 
differential cross section parameters was itself parameter- 
ized and then fed back into the final fit. The fit was then 
re- minimized. Free parameters and constraint terms were 
included to allow the energies to vary somewhat in the 
fit. This approach did not allow the fit to explicitly ex- 
plore the effects of these perturbations on subtleties such 
as changes in center-of-mass angles of detectors. We em- 
phasize "explicitly" because such subtleties were in fact 
accounted for when originally parameterizing the effects 
of the perturbations. For this reason, all freedom in the 
final results arising from such effects was still quantified 
and reflected in the error matrix. The increase in the 
uncertainties calculated by this fitting procedure was in 
excellent agreement with quadrature addition of the tab- 
ulated energy-related uncertainties to the original uncer- 
tainties in the fit. 



4- Peak-Fitting Uncertainties 

Some ambiguity existed in the fits of the Rutherford 
backscattering spectra for deuterons scattering from gold 
at the lowest energies. It was generally unnecessary 
to account for correlations arising from uncertainties in 
monitor detector yields. However, because of the sys- 
tematic nature of this normalization peak and the non- 
negligible uncertainty in the fit, this rationale became 
invalid. We added extra normalization parameters for 
these data in the global fit. Corresponding constraint 
terms were added to the x 2 sum constraining these pa- 
rameters to be near one. The uncertainties assigned to 
the constraint terms were 3% for the 120 keV normal- 
ization parameter and to 1% for all other energies at or 
below 320 keV. 



IV. RESULTS, IMPLICATIONS, AND 
CONCLUSIONS 

Here we present our data and address the effect these 
data will have on Big-Bang Nucleosynthesis calculations. 
The cross-section results obtained from the parameteri- 
zation described in Eq. are given in Table II I II The 
coefficients of the Legendre expansions of the differential 
cross-sections are given in Table IIVI The errors shown 
are la equivalents which include statistical and back- 
ground uncertainties, energy uncertainties discussed in 
Sec. IIII E~3l fitting uncertainties described in Sec. IIII E4l 
and which account for scatter in the data sets used to de- 
rive the total cross sections as explained in Sec. IIII E 21 
The uncertainties are generally around the 2% level ex- 



cept at the lowest energies where peak fitting and en- 
ergy uncertainties become somewhat larger. An overall 
scale error of 1% is estimated arising from uncertainties 
in the p-d elastic-scattering cross section to which these 
data were normalized. This normalization error is not 
included in the uncertainties. Finally, the error matrix 
for the total cross-section parameters is given in Table IVl 
The complete error matrix for the Legendre coefficients 
is too large to present here and is not relevant to our 
primary goal of providing BBN inputs. 

The cross sections of Table IIIII are plotted in Fig. [7| 
along with the continuous parameterization mentioned 
in Sec. IIII Dl and with the recent cross section compila- 
tion of Cyburt 26]. Some scatter within the uncertain- 
ties is detectable, which is expected since our statisti- 
cal and systematic errors are of comparable magnitudes. 
The continuous parameterization appears to be system- 
atically somewhat higher than the results of the discrete 
fit. This is not surprising since cross sections for the 
lowest four energies are all normalized directly to the 
differential cross sections at 480 keV. These data thus 
inherit all of the uncertainty of the 480 keV points, but 
in a systematic manner. The continuous fit is then free 
to systematically renormalize these points within con- 
straints. This systematic uncertainty is quantified in the 
error matrix of the discrete fit. 

The Cyburt compilation agrees well with the present 
data except at the high-energy end the 2 H(d,p) 3 H reac- 
tion cross sections. At the lowest energy we can com- 
pare directly to the high-precision data of Brown and 
Jarmie [23, which falls 11% and 8% below the present 
work for 2 H(d,p) 3 H and 2 H(d, n) 3 He respectively, corre- 
sponding to discrepancies of 2.1 a and 1.6 a respectively 
when systematic and statistical uncertainties of both ex- 
periments are included. The data of Greife et al. [28| . 
having uncertainties of only about 2.8%, agree well with 
our results at higher energies, up to 256 keV, for both 
reactions. 

There have been several published compilations and 
analyses of BBN data and network calculations and we 
cannot review or even acknowledge all of the major efforts 
here. Most, including Refs. 6, 29), |3(J, |3lJ , have focused a 
significant amount of attention on consistency of cosmo- 
logical observables including primordial abundances and 
CMBR observations. Many, including some which have 
now acquired benchmark status, have used Monte Carlo 
techniques to analyze the relationship between the uncer- 
tainties in the reaction data to the uncertainties in the 
astrophysical constraints [8l l32L l33| . The work of the lat- 
ter produced estimates of energy regions having the high- 
est sensitivities to cross-section data, estimates which in- 
spired the present work. Finally there are even some 
works that have used traditional first-order error propa- 
gation to understand the effects of the data on network 
predictions. For example, Fiorentini et al. |34| achieved 
impressively good agreement with the computationally 
expensive techniques used elsewhere and have the ad- 
vantage of producing simple functional understandings 
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of various uncertainty relationships via tabulated deriva- 
tives. 

The NACRE compilation has become widely ac- 
cepted as a standard in adopted values for reaction- 
rates of nuclear reactions of astrophysical significance. 
NACRE does not address network calculations since it 
is intended as a broad resource for general astrophysical 
use. It has though become the data source of choice for 
several BBN analyses pfl IM Isfij . 

The NACRE reaction-rate values cannot be compared 
directly to the cross-section data of the present work. 
The reaction rate is an integral of the total cross section 
weighted by a Maxwellian distribution and taken over all 
energies. In order to calculate the reaction-rate integral 
for the purpose of this comparison, we used the cross sec- 
tion results from the continuous parameterization of the 
present work. For energies given in MeV and cross sec- 
tions in barns, the reaction rate N a (av) in cm 3 mol _1 s _1 
is given by 

/>oo 

N a {av) = 3.731 xl0 10 /i" 1/2 ^ 3/2 / oE eT 11 ' 604 * dE , 

Jo 

where [i is the reduced mass in amu and Tg is the tem- 
perature in units of 10 9 K. 

For the integral to converge at temperatures relevant 
to BBN, T 9 < 2, it can be cut off at E c , m . of about 2.5 
MeV but not much lower. Thus to calculate reaction 
rates accurately, we must include data at energies higher 
than those of the present work, although data at these 
energies contribute relatively little to the integral at the 
relevant temperatures. For this purpose and to make 
a fair judgment of how the new data compare to the 
NACRE compilation, we have used the same data at high 
energies as those used in NACRE compilation, the data 
of Schulte et al. ^S%. We have included this data up to 
£ c m = 2.75 MeV. 

These data were included in the continuous parame- 
terization by using the total cross sections to appropri- 
ately constrain the zero order Legendre coefficients at the 
high energies. The data of Ref. [37| have scatter which 
is significantly larger than their quoted statistical errors. 
To prevent this data set from unduly constraining the 
fit, we have multiplied their statistical uncertainties by 
5, leaving the smallest uncertainties in these data still 
below one percent of the value. This multiplication is 
a coarse procedure but suffices for the present purpose. 
The results of these fits are shown in Figs. El and The 
2 H(<i,n) 3 He reaction curve does not follow the Schulte 
data impressively well. This is probably due to a limi- 
tation of the parameterization in representing changes in 
curvature over these large energy regions. It is not a sig- 
nificant problem for the purpose of constraining the tail 
of the reaction-rate integral. 

The NACRE compilation gives coefficients for a 
quadratic cross-section function which was used in calcu- 
lating the rates for the 2 H(d, n) 3 He and 2 H(d, p) 3 H reac- 
tions at " low energies" . NACRE does not specify exactly 



what is meant by low energies. We have used these cross- 
section functions for energies below the present work in 
order to complete the integrals. As was the case for the 
high energy data, this low energy data contributes only a 
small amount to the integral for temperatures significant 
to BBN. 

Using these total cross-section curves, we then calcu- 
lated the reaction rates by computing the right hand 
side of Eq. |2] with a high-energy cutoff on the inte- 
gral at -E c .m. = 2.5MeV. The results are compared to 
the NACRE reaction rates in Figs. an d El Judging 
by plots of abundances as a function of temperature for 
network calculations of Nollctt and Buries [g] , all signifi- 
cant standard BBN occurs at temperatures near Tg = 1. 
At this temperature the rates derived using the present 
work are about 7% higher than the NACRE rates for 
both 2 H(d, n) 3 He and 2 H(rf,p) 3 H . These discrepancies 
are systematic, remaining at similar levels for a large 
range of temperatures. Our derived 2 H(d,p) 3 H rate is 
well beyond the NACRE quoted upper bounds, and for 
2 H(d, n) 3 He our value is marginally within the NACRE 
bounds. 

It is important to emphasize that there is nothing in- 
herently controversial about disagreement between the 
present work and NACRE or any other compilation. The 
NACRE curves are derived from a limited set of data 
and in fact include no data between Ed = 325 keV and 
Ed = 2.0 MeV. Furthermore, the data of Krauss et 
al.. were the only data used in the NACRE compila- 
tion in the energy range of 250 to 325 keV, and although 
these are respectable data, they have total uncertainties 
of about 8% (see Fig. [I}. The NACRE compilation incor- 
porates all data in these plots with the exception of the 
Ganeev data. Our good agreement with Cyburt [2(| is 
likely due to his inclusion of the Ganeev data set which, 
although it has large uncertainties, also agrees well with 
our data and covers the energy range with the least data 
available. There was indeed a previous shortage of data 
and our results will greatly contribute to the cross-section 
information available in this energy region. 

Many recent compilations such as that of Dcscouvc- 
mont et al. 39] used in the analysis of Coc et al. l36|. as 
well as the compilation described by Serpico et al. [40(, 
have a focus on BBN applications and strive to improve 
on the extrapolation techniques of the NACRE compila- 
tion in the BBN energy regions. However, these compila- 
tions suffer from a lack of data in the same energy region 
as NACRE. In spite of this, Serpico et al. quote unusu- 
ally low uncertainties at the levels of 1.3% and 1% for 
the 2 H(d,n) 3 He and 2 H(d, p) 3 H reaction rates respec- 
tively. They do however acknowledge the severe lack of 
data and " strongly recommend a new experimental cam- 
paign" which we have now furthered. 

Us ing the logarithmic derivatives tabulated in 
Ref. |3J], we can estimate the differences in primordial 
abundances calculated using our new data vs. using the 
NACRE compilation. These derivatives of abundances 
with respect to reaction rates are given in tables of coef- 
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ficients of polynomial expansions in the baryon to pho- 
ton ratio 77. For our estimates we use a value of 77 of 
6.14 x 10 10 from Ref. 6] which was derived directly from 
an analysis of recent WMAP data performed in Ref. • 
The values of the logarithmic derivatives calculated for 
77 = 6.14 x 10 10 are given in Table EH 

It is immediately evident that a roughly 7% change in 
the cross sections will have essentially no effect on the 
primordial abundance of 3 He . The relative signs of the 
derivatives for the two reactions have important signifi- 
cance. Although the derivatives shown for the 3 He abun- 
dance are of significant magnitude, the changes in the 
abundance induced by the changes in the two reactions 
nearly cancel and the net effect is small. However, be- 
cause the two D /H derivatives have the same sign as do 
the reaction rate changes inferred from the present work, 
the total change estimated in the deuterium abundance 
is roughly —7%. The change calculated for 7 Li/H is +5% 
and is almost entirely attributable to the changes in the 
2 H(d, n) 3 He reaction rates. 

The lithium abundance observations are currently in 
severe disagreement with the Standard BBN results [2|]. 
This 5% increase in the predicted lithium abundance 
makes an unresolved problem slightly worse. Currently 
the 7 Li abundances predicted using CMBR values of 77 
are at least a factor of 2 higher than the observed abun- 
dances 29] . The 5% change is small compared with the 
uncertainties plaguing this comparison but it certainly 
does not improve the hopes of finding agreement between 
standard BBN and observed 7 Li abundances. We look at 
this not as a failure but as an opportunity to explore pos- 
sibilities of non-standard Big Bang Models and thus new 
physics. 

The current observational value of the primordial deu- 
terium abundance, D/H = 2.78±g;|| x 1CT 5 0, ob- 
tained from QSO absorption spectra, is in good agree- 
ment with predictions of BBN using 77 from the CMBR 
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data, D/H = 2.561^ x 10" 5 |29| or D/H = 2.55 
10 -5 However, as more data have arrived, the sta- 
tistical uncertainty estimates on the abundance obser- 
vations have not held up; the uncertainty is now dom- 
inated by scatter in the data. The change which we 
predict in the D/H value derived from BBN+CMBR 
will strain the current agreement slightly, but the results 
should still be well within the current mutual uncertain- 
ties. The reduced uncertainties of the 2 H(<i,n) 3 He and 
2 H(d,p) 3 H cross-sections from the present work, future 
improvements in precision of the 2 H(p,7) cross-section 
measurements, and more QSO observations may soon 
provide significantly more stringent tests of this compar- 
ison, again yielding insight into the validity of the details 
of the standard BBN model. 

Although other aspects of the BBN picture are still 
limiting the comparisons of theory and experiment to 
around the level of 10% or worse, the new data of the 
present work will pave the way for future developments 
in precision cosmology. As other measurements are im- 
proved, BBN predictions will become sensitive to many 



details of the model including possible inhomogeneities 
and neutrino properties. The new data should put to 
rest concerns and claims that the d-d reaction rates may 
be inaccurate, and should ultimately result in significant 
modifications to the reaction rates used for current astro- 
physical applications including BBN. The present data 
verify the limited, not always trusted, and often over- 
looked data sets which previously existed in this region of 
energy ^1 an d with roughly 2% to 3% uncertainties, this 
work represents a significant improvement in precision, 
which along with improvements in other cross-sections, 
will translate directly into reductions in uncertainties of 
BBN predictions. These data will truly help to usher in 
the " new era" |(| of precision cosmology. 
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TABLE III: Integrated cross-section results for discrete fit. 
The corresponding S-factor error matrix is given in Table [V] 
The la uncertainties include statistical uncertainties from all 
parts of the present work as well as fitting uncertainties dis- 
cussed in Sec. 1111 El 

E d (keV) a-2 H(d>p) 3 H (mb) q-2 H(dj „ ) 3 Hc (mb) 

646.1 74.55 ± 1.07 84.34 ± 1.59 
557.3 68.02 ± 0.96 80.22 ± 1.18 

470.2 58.47 ± 0.80 69.55 ± 1.04 
379.2 53.09 ± 0.94 62.60 ± 1.12 
314.6 45.77 ± 1.06 56.01 ± 1.33 
232.8 37.65 ± 0.82 43.13 ± 0.96 
172.8 29.55 ± 0.88 33.03 ± 1.01 
112.2 19.77 ± 0.85 21.07 ± 0.92 
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FIG. 7: (Color online) Integrated cross-section results for 2 H(d,p) 3 H and 2 H(rf, n) 3 He . The points are the results from the 
discrete parameterization of Eg. 1211 and the curves are the results of the continuous parameterization. The error bars are la 
uncertainties including statistical uncertainties as well as fitting and energy uncertainties. Details are discussed in Sec. 1111 El 



TABLE IV: Legendre coefficients for discrete fit of 2 H(d,p) 3 H and 2 H(d, n) 3 He data. The la uncertainties include statistical 
uncertainties from all parts of the present work as well as fitting and energy uncertainties discussed in Sec. ITTTEl The coefficients 
are normalized to the Po coefficients. 



£ d (keV) 



646.1 
557.3 
470.2 
379.2 
314.6 
232.8 
172.8 
112.2 



l>2 



In 



0.786 ± 0.004 0.192 ± 0.005 
0.737 ± 0.004 0.199 ± 0.006 
0.665 ± 0.007 0.075 ± 0.010 
0.629 ± 0.007 0.094 ± 0.008 
0.536 ± 0.005 0.079 ± 0.010 
0.472 ± 0.005 0.013 ± 0.010 
0.463 ± 0.006 0.046 ± 0.007 
0.358 ± 0.012 0.024 ± 0.008 



"H(d,n) 3 He 



l>2 



In 



0.865 ± 0.032 0.011 ± 0.030 

0.856 ± 0.011 0.156 ± 0.013 

0.798 ± 0.018 0.048 ± 0.025 

0.830 ± 0.016 0.112 ± 0.019 

0.721 ± 0.008 0.029 ± 0.016 

0.663 ± 0.007 0.011 ± 0.011 

0.697 ± 0.017 0.059 ± 0.022 

0.578 ± 0.010 0.016 ± 0.008 
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TABLE V: Error matrix for discrete fit. the elements shown are for the zero order Legendre coefficients, the ao's. The two 
reactions 2 H(d,p) 3 H and 2 H(d, n) 3 He are labeled by p and n respectively and by their lab energies in keV. 
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FIG. 8: (Color) Integrated cross-section result for 2 H(ci,p) 3 H 
for fit of continuous parameterization including integrated 
cross-section data of Schulte et al. [37l| The statistical un- 
certainties of Ref. [37| were multiplied by 5 in the figure and 
the fit. 
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FIG. 9: (Color) Integrated cross-section result for 
2 H(rf, n) 3 He for fit of continuous parameterization including 
integrated cross-section data of Schulte et al. [37l | The sta- 
tistical uncertainties of Ref. [3^1 were multiplied by 5 in the 
figure and the fit. 
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TABLE VI: Logarithmic derivatives of abundances Y with 
respect to reaction rates R calculated from Ref. |34l . 



Abundance Ratio 
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FIG. 10: (Color online) Comparison of 2 H(rf,p) 3 H reaction 
rates of present work with the NACRE rate compilation. The 
results of the present work use the low energy cross-section fit 
of NACRE in the rate integral for E c . m . < 0.05 MeV. Data 
of Ref. |37| are used to constrain the cross-section fit of the 
present work at energies above -Ec.m. = 0.5 MeV. 
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FIG. 11: (Color online) Comparison of 2 H(ii, n) 3 He reaction 
rates of present work with NACRE rate compilation. The 
results of the present work use the low energy cross-section 
fit of NACRE in the rate integral for E c . m . < 0.05 MeV. Data 
of Ref. J37J are used to constrain the cross-section fit of the 
present work at energies above E c . m . — 0.5 MeV. 



